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Abstract 

Microscopic variations in composition or structure can lead to nanoscale inhomogeneity in super- 
conducting properties such as the magnetic penetration depth, but measurements of these prop- 
erties are usually made on longer length scales. We solve a generalized London equation with a 
non-uniform penetration depth, A(r), obtaining an approximate solution for the disorder-averaged 
Meissner screening. We find that the effective penetration depth is different from the average pen- 
etration depth and is sensitive to the details of the disorder. These results indicate the need for 
caution when interpreting measurements of the penetration depth and its temperature dependence 
in systems which may be inhomogeneous. 



PACS numbers: 74.62.En,74.20.-z,74.62.Dh,74.81.-g 



I. INTRODUCTION 



The penetration depth and its temperature dependence are important characteristics of 
any superconductor and are considered key to determining the momentum space structure of 
the order parameter.- - - The possibility of disorder in exotic superconductors is well known, 
but analyses performed to date have concentrated on the effect of disorder-induced scattering 
on the momentum space structure of the gap." This paper is motivated by the possibility 
that disorder may lead to nanoscale real space variation and the associated need to model 
the relationship between such spatial variation and properties that are measured on longer 
length scales. We address how inhomogeneity in the penetration depth may affect bulk 
measurements of the penetration depth for methods that rely on Meissner screening and can 
be analyzed by solutions to London's equation. In particular, we show that the measured 
result is not simply given by the average value of the penetration depth, but is affected by 
the statistical structure of the spatial variations in the penetration depth. 

Many superconductors are created by chemical doping of a non-superconducting parent 
compound. In these systems the inherent randomness of the doping process may give rise 
to an inhomogeneous superconducting state. The importance of this effect will be deter- 
mined by the characteristic length over which the dopant atoms affect the superconductivity. 
Even in the most ordered material, there will be binomial fluctuations in the total number of 
dopants in a given region. In general, one does not expect significant spatial variation in ma- 
terials that are weakly correlated and can be described by a rigid band model. For example, 
disorder is largely irrelevant in classic metallic superconductors, due to their long coherence 
lengths and weakly correlated nature.- In contrast, the cuprates are doped insulators with 
a coherence length on the scale of the lattice. They are known to have nanoscale disorder 
in their superconducting properties, as seen by scanning tunneling microscopy. 9 Similar gap 
maps have been observed in the iron pnictide family^— and in disordered titanium nitride 
films close to the superconductor to insulator transition.— 1 ^ 

Materials with intrinsic disorder present two separate challenges. Understanding how the 
random doping process gives rise to local superconducting properties, such as the penetration 
depth or local density of states, requires a microscopic model. But even with such a model, 
we still need to make the connection between the local superconducting properties and bulk 
measurements. The manner in which local superconducting properties relate to the observed 
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properties will differ from experiment to experiment. For instance, a measurement of the 
heat capacity will return the total heat capacity of the macroscopic sample, so the inferred 
specific heat capacity will be a volume average over the sample. In contrast, we might expect 
the thermal conductivity response to be dominated by a percolation path connecting regions 
with small local gap, A(r), or large local density of states. 

Here we focus on the penetration depth, A, as measured by screening of the magnetic 
field, including both resonant cavity frequency shift measurements at radio frequencies^ and 
the local probes of Magnetic Force Microscopy^ and Scanning SQUID Susceptometry.— 
These methods measure A(T) by detecting the response magnetic field generated by the 
superconductor due to an applied field, and can be analyzed using the London equation. 
Thus, we can model the effect of inhomogeneity by solving the London equation with A(r) as 
a random function of position r. Our goal is to find a new equation for the disorder-averaged 
magnetic field, as this will determine the measured response. Here, we work in the limit of 
small fluctuations to find an approximate equation for the disorder-averaged magnetic field, 
as this will determine the measured response. 

II. STOCHASTIC LONDON EQUATION 

To understand the measured penetration depth when A(r) is a random function of po- 
sition, we calculate the disorder-averaged magnetic field response to obtain an effective 
penetration depth. For isotropic and local superconductors in three dimensions, the static 
magnetic field h(r) is given by the London equation with A(r) a function of position. The 
correct form^ of the London equation when the penetration depth is non-uniform is: 



which is derived from the second Ginzburg-Landau equation in the London limit.— We 
parametrize the penetration depth as the average value plus a fluctuating term: 




(1) 




(2) 



so that (A(r)) = A. Then Eq. CD becomes: 



(L + M l + M 2 ) h = 0, 



(3) 
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where: 



M 2 



Mi 



L 



1 - A 2 V 2 J, 

-2A 2 £V 2 I + 2A 2 V£ X V X , and 
-A 2 £ 2 V 2 J + A 2 V£ 2 x Vx. 



(4) 



Here I is the identity tensor, and the "dangling curl" is understood to operate on a vector to 
its right. The terms are grouped so that Mi is first-order in £, M 2 is second-order in £, and 
L gives the unperturbed London equation. We will work in the limit of small fluctuations, 
£(r) <C 1, so that Mi + M 2 is a perturbative term in Eq. [31 

Our method of solution comes from the similarity of the Helmholtz and London equations. 
The Helmholtz equation, which governs wave propagation, becomes the London equation 
when the wavevector is purely imaginary. Thus our problem is related to the propagation of 
waves in a random medium, and we can build upon a large and multidisciplinary literature 
devoted to this challenge.— 1 ^ The paper by Karal and Keller— is particularly relevant, 
because it retains the vectorial nature of the problem, rather than simplifying to a scalar 
wave equation. 

We now derive, from Eq. [31 an approximate equation for the disorder-averaged field (h). 
Applying the inverse of L to both sides: 




(5) 



where L h = 0. Solving for h: 



h= [1 + L- 1 (M 1 + M 2 )] h ( 



(6) 



assuming the inverse exists. Averaging both sides: 




(7) 



where ho comes outside of the average because it is non-random. Solving for h : 




(8) 



Since we assume small fluctuations, we can expand the term inside the average: 
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Averaging and expanding again: 

(1 - L- 1 (M 1 L~ 1 M 1 ) + L-\M 2 )) (h) = h , (10) 
since (Mi) = due to Eq. [2J We then apply L to both sides: 

(L - (M 1 L~ 1 M 1 } + (M 2 )) (h) = 0, (11) 
which yields the average field to second order in £. 

III. RESULTS 

We first evaluate the averages in Eq. [TTJ giving us an equation for (h) in terms of the 
penetration depth correlation function, (A(r)A(r')). We then consider two specific cases for 
the correlation function and numerically evaluate the effective penetration depth for a range 
of parameters. 

A. Evaluating the Averages 

We will solve Eq. [TT]for a single Fourier mode of (h(r)) = h e lk r , then derive an equation 
for k that yields exponentially decaying solutions consistent with Meissner screening. 
First, we evaluate (M 2 ): 

(M 2 ) = -A 2 (£(r) 2 ) V 2 / + A 2 ( V£(r) 2 ) x V X . (12) 

We introduce the correlation function R(r,r') = (£{r)£(r')), which is a function only of 
\r— r'\ if £ is stationary and isotropic. Then we see that (£ 2 ) = -R(O) and (V£ 2 ) = V(£ 2 ) = 0, 

so 

(M 2 )(h) = A 2 k 2 R(0)he lkr . (13) 

We now evaluate the remaining average, <(M 1 L _1 M 1 ), in three stages to derive Eq. [2TJ 
First we expand the differential operations, then evaluate the disorder average. The last 
stage is to evaluate the integral. We will then combine this integral with Eq. [T2] to solve 
Eq. HJ 

The average to evaluate has the form: 

(MiL _1 Mi)(/i) = /"dr' / ^Mi(r)G(r-r')Mi(r' / ))(/i(0). (14) 
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The Green's function is the solution to (1 — A 2 V 2 ) G(r, r') = 5(r — r'), and is: 

1 1 



G(z) 



-z/A 



(15) 



A 2 4irz 

Here, we have defined z = r — r' . 

We now expand the differential operations in Eq. HU We do this in two segments, first 
with derivatives at r, then with derivatives at r'. The first is: 

Mx(r)G(r - r')v{r') = [-2A 2 ^(r)V^ + 2A 2 V r £(r) x V r x]G(r - r>(r') 

= 2£v{r') [5{z) - G{z)] + 2A 2 V£ x [V r G(r - r')} x v{r'). (16) 

The second, which was represented by v(r') above, is: 

w(r') = Mx{r')(h(r')) = e ikr ' \2A 2 k 2 £(r')h + 2tA 2 VA(r') xkxh 

= 2A 2 e lkr ' | [k 2 i(r') - i V P /f (r') -k]h + i [V r /f (r') • h] fc}. (17) 
Combining Eqns. [TO] and [T71 we obtain: 

Mx{r)G{r - r')Mx{r'){h{r')) = 2AV fcV jf (O [5(2) - G(,s)]2A; 2 h 

+ £,{r)V r >£(h ® k - k ® h)2i[5(z) - G(; 



-2A 2 G(^)(A- 1 + z' 1 ) 



f (r')V P f (/i ® z - z ® h)k 2 



+ i(h®z-z®h) (V r £ ® V r >£)k + i(z®k-k®z) (V r £ ® V P 'f)/i 



(18) 



where we use ® to indicate the tensor product. 

To perform the disorder average in the second stage, we need various derivatives of the 
correlation function R(z): 

(C{r) VA) = V r / R(\r - r'Q = -z R{z), 
(C(0 V r ^) = V r R(\r -r'\)=z R(z), and 

(V r £V r ,0 = VV'R(\r-r'\) = 

where the overdot indicates differentiation with respect to z, and / is the identity tensor. 
Then averaging Eq. [TBI gives: 

(Mx{r)G{r - r')Mx{r f )){h{r')) = 2AV fcr '{ [A(z) + 2k 2 R(z)5(z)]h 

- [B(z) + 2iR(z)6{z)](k®h-h®k)z-C{z)h(z®z)\, (19) 



R T , . R\ 
—I+z®z[R 

z \ z J 



with the scalars A, B, and C given by: 



A(z) = 2k 2 R(z)A 2 G(z)(A- 1 + z' 1 ) - R{z) G(i 
B{z) = 2i \r{z)A 2 G{z) (A^ 1 + z- 1 ) - R{z) G{z) 
C(z) = 2A 2 k 2 R(z) G(z) (A' 1 + z' 1 ) . 



and 



The final stage in evaluating Eq. [HI is to perform the integral over r' . We first change 
variables from r' to z, then integrate over the orientation of z. Using the relations 

sin(kz) 



dz e- lk z = A-k- 



kz 



F(k,z) 



dz ze ik z = k - d k F, and 



(20) 



dz z ® z e 



-1 



d,F 



I + k®k[d 2 k F - 



d k F 



we find that Eq. [TH evaluates to: 



J rfr , ^M 1 (r)G(r-r / )M 1 (r , ))(^(r / )) = 4A 2 k 2 e lk r {[X + R(0)} ^ + (21) 



The functions X and Y are given by: 



X(k) = J dzG{z){ A 2 (A- 1 +z~ 1 ) R(z 2 F+k- 1 d k F)-Rzk- 1 d k F +Rzk~ 1 d k F-Rz 2 F \, 

"CO 



Y(k)= I dz G(z){A 2 (A- 1 + z- 1 ) Rzk- 1 d k F + R(d 2 k F-k- 1 d k F) 



- Rzk~ x d k F 



(22) 



We require the average magnetic field to have V-/i = 0, which means that k ■ h = 0. We 
now collect our results from Eqns. [21] and [131 an d insert them into Eq. [TT1 



he 



ikr 



1 + A 2 k 2 (l-3R{0)-AX) 



0. 



(23) 



We are interested in solutions consistent with Meissner screening, so we require that k be 
positive and purely imaginary. Then the field decays on a length scale A e g = |, which we 
identify as the experimentally measured penetration depth. To calculate A e e, we will solve 
the equation: 

= 1 - 3R(0) - AX. (24) 

A 2 



Inserting k = ^- into our equation for X, we get: 

X = 4vr ^ dz G(^)sinh (f) | A 2 (A-^" 1 ) 

+ 4vr^ dz G(,2)cosh^-^-^ |a 2 (A _1 + z~ l ) 

Valid solutions for X c g will require the integral for X to converge and Eq. [2H to have 
solutions. 

B. Correlation Function 

A full solution of the disorder-averaged magnetic field, (h), requires knowledge of the cor- 
relation function R(z) and hence requires not only a detailed knowledge of the composition, 
structure, and disorder of the sample, but also a microscopic model to locally determine 
the superconducting properties from that structure. Without guidance from microscopic 
calculations, we will use the Matern one-parameter family of correlation functions^ 2 to tune 
the smoothness, as well as the magnitude and correlation length, of the penetration depth 
fluctuations. Handcock and Wallis™ parametrize the Matern class of covariance functions 
as: 




where K v is a modified Bessel function of the second kind and T(z) is the Gamma function. 
The intercept at zero separation is the normalized variance of the penetration depth, -R(O) = 
a \/ W 2 = ((^ 2 ) ~~ (^) 2 )/(^) 2 5 an d quantifies the magnitude of the inhomogeneity in A(r). 
The correlation length, /, controls the size of the fluctuations in A(r). The parameter v 
controls the smoothness of A(r). Larger v gives a smoother random field, since it is \v\ — 1 
times mean squared differentiable, where [•] is the ceiling function.— 

Two members of the family deserve specific mention. When v = 1/2, Eq. [261 reduces 
to the exponential correlation function, R(z) = R(0) exp(— zy/2/l^j , which is the correlation 
function of a Markov process in one dimension. The integrals for X in Eq. |25] diverge 
when R(z) oc e~ z , making the case v — 1/2 invalid. In the limit where v — > oo, R{z) — > 
R(0) exp(— z 2 /l 2 ), which is labeled the squared exponential correlation function, to prevent 
confusion with the Gaussian probability distribution. This correlation function gives the 
smoothest possible A(r) that can be described within the Matern covariance family. 
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-RX 3 



off 



+RX 3 eS -RzX 



cff 



-.ra^ + i^A 2 



cff 



-RzX 2 



eif 



(25) 



C. Squared Exponential Correlations 



We now consider the case of squared exponential correlations, R(z) = R(0) e~ z2 ^ 2 . In 
Fig. [T] we plot four realizations of a normally distributed penetration depth with squared 
exponential correlations, illustrating the effect of the two parameters / and R(0) on A(r). 
Evaluating Eq. [251 gives: 

A- 



X = -R(0) / dz e~ z/A e~ z /l sinh 
Jo 

POO 





A e ff 


Geff) 


A 2 " . 



(l + 2 £)+2z-j (1 + 2^) 



+ 2R(0) / dz e~ z/A e~ z /21 cosh 
Jo 

All of these integrals converge, so we evaluate X as: 
X = R(0) 
+ R(0) 



\ 2 

A eff 

A 2 



: - (1 + 2^) +2z : 



1 



A 

¥ 



(27) 




2l 2 A\ e{i + AA 2 \ 2 eS exp 



4PA 2 

- (V + 2/ 2 AA cff + 4A 2 A c 2 ff ) exp 



1 



1 
A 

1 

Aofr 



1 

A e ff 



erfc 



erfc 
2 V A 



1 1 

-r + t — 

A Aeff 



1 

Aeff 



(28) 



After inserting Eq. into Eq. [2H we solve for A e s over three decades in the corre- 
lation length, /, and in the disorder variance, -R(O) (Fig. [5]). At large correlation length 
the effective penetration depth is larger than the average value, representing suppressed 
Meissner screening. Conversely, at small correlation length the effective penetration depth 
is smaller than the average, indicating enhanced screening. The separatrix, where A c ff = A 
for all values of R(0), occurs near I = 1.643A. Note that the system is not symmetric 
about the separatrix, although it becomes more symmetric as R(0) —> 1. This is true 
for both linear and logarithmic spacing around the separatrix. In other words, neither 
|A eff (/ s + A/) - (A) | = | A c ff (/ s - A/) - (A) | nor |A eff (a/ s ) - (A}| = |A cff (Z s /a) - (A)| are true, 
where l s denotes the separatrix, and a is an arbitrary positive real number. As expected, 
A e ff — > A as R(0) — > 0. Yet even at small disorder, A e ff has variations on the one percent 
scale, shown by the contours in Fig. |2J As we will discuss below, sub-percent variations of 
A e ff could be significant in the context of a typical measurement of AA(T). 

The trends in A e ff can also be seen in Fig. [21 where we plot A e ff/A vs. R(0) at fixed 
correlation length. All three curves taper to A e s = A as the magnitude of disorder decreases. 
At large correlation length, in this case I = 10A, A c g increases by ten percent when R(0) = 



0.02. The effect at small correlation is more modest, but still reaches nearly ten percent by 
the time R(0) = 0.1 when / = 0.1A. 

The penetration depth has a temperature dependence that it inherits from the underlying 
disordered superconducting state. It is natural to expect that R(0) and I will have a tem- 
perature dependence of their own, which will create a temperature-induced change in \ cS . 
This change contributes to any measurement of A(T), but is not related to the gap structure 
in momentum space, because it arises from the spatial arrangement of the superconducting 
state. If we neglected the spatial variation of A we would erroneously attribute the entire 
temperature dependence to the order parameter. 



D. General Matern Correlations 



To understand the impact of the smoothness of A(r) on the measured penetration depth, 
A e ff, we now consider the general case of Matern covariance. Recall that the parameter v 
controls the smoothness of the penetration depth. With the correlation function defined by 
Eq. [261 we evaluate Eq. 



X 



R(0) 



2^- 1 rM 



+ Aul 2 (A 2 + Az + X 2 cS ) (2y/Vj) i^_i(2V^ 



u-l 



+ 16z/ 2 A 2 ff A(A + z) (2v^yj K u . 2 [2y/D- 



u-2 



+ 



R(0) 



dz e z/A cosh 



4^A c 2 ff 
A 2 / 4 



/ 2 : ( 2^)^^-1(2^ 



+ 4uA(A + z)z (2^1)" 2 K u _ 2 (2y/V 



■ (29) 



These integrals can be evaluated using equation 6.621.3 in Gradshteyn and Ryzhik: 24 

POO 

/ x"- 1 e- ax K v (/3x)dx = 
Jo 



(2/3)^ Y^ + u)Y^-u) 



1 1 oc — (3 

2 Fi \ /i + v, v + 5 ; fj, + |; 



a + (3 



(30) 



which requires Re/i > \Kei/\ and Re (a + 0) > 0. The function 2 Fi(a,b; c; z) is Gauss' 
hypergeometric function. Using the integral in Eq. [30] to evaluate Eq. [291 we find the 
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constraints 



3 A c g I 

v > - and — — > 7 — 

2 A I + 2Ay/u 



(31) 



The full solution for X is then: 

R(0)V^ (jvY Ace f A 2 ff + A 2 r(2z/-1) 

i» V^ 2 / a [ 2A i> + ±) 
r(2i/) 



-(2u 



+ 



r(*+§; 



A e ff + A _ 2j , 



2A 



a- 2u F2(D) 



A fiff - A 



2A 



b- 2v F2(+) 



+ 



+ 



i r(2z/ + 1) 
a i> + §) 

A c 2 ff Ar(2z/-3) 
4 T(u-l) 
A eff r(2z/-2) 



-(2j/+1) 



+ 



-(2i/ 



F3(D) - r (2 " +1) F3(*) 

- 3) F4(n) - r (2 "- 3) F4(*) 



4 r(i/ + |) 



(A eff + A)a-^- 2 )F5(D) - (A eff - A)&-< 2 "- 2 >F5(*) 



+ 



A eff r(2^-i) 
2 r(i/ + |) 



,-(2f- 



^F6(D) + r (2l - 1} F6(*) 



(32) 



where we have introduced the variables 

1 



and functions 



A^A cfr ^ I ' 



□ 



1 

A 



— + ^. 

A e ff / 



/(A c ff + A) - 2A c ffAy^ 

/(A eff + A) + 2A eff A07' 

/(A - A c ff) + 2A c ffAy^ 

/(A - Aeff) - 2A eff A x /i7' 





= 2 F 1 (2v-l,v-\- 


z/ + 


!;•) 


F2(-) 


= 2 F l (2v,is-±;v + 


2' / 




F3(-) 


= 2^(21/ + + |; 


z/ + 


§;•) 


FA(-) 


= 2^1(21/- 3, i/-§; 


z/ — 


§;•) 


Fh{-) 


= 2^(21/ -2, i/- |; 


z/ + 


I;-) 


F6(-) 


= 2 Fi(2i/ -l,i/-f; 


z/ + 


§;•) 



and 



Inserting this expression for X into Eq. [24J, we can solve for A e ff after choosing a value 
for the smoothness parameter v. In Fig. HJ we have chosen v — 2, close to the lower bound 
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of I required for convergence of X. The results are almost identical to the case of squared 
exponential correlations (Fig. [2]); evidently A e fj is not much affected by changes in the 
smoothness of A(r) for the Matern family of correlation functions. The qualitative features 
of interest to us are still present: there are regions of enhanced screening and regions of 
suppressed screening, the effect grows on increasing the variance of A(r), and changes in A e fr 
at the one percent level persist down to small disorder. Quantitatively, the results in Figures 
H] and [2] differ by five percent in the region near I = 1 and -R(O) = 1, where the difference is 
largest. 

IV. DISCUSSION 

The measured A(T) in a non-uniform superconductor will be determined by both the 
momentum space gap structure and the real space variations of the penetration depth. We 
calculate the influence of spatial fluctuations in the penetration depth by solving the stochas- 
tic London equation in the limit of small fluctuations. This gives an equation (Eq. 124]) for 
the disorder-averaged magnetic field in terms of the penetration depth correlation function. 
We then solve this equation for two example correlation functions to find A e fr, the decay 
length of the disorder-averaged field, which we identify as the penetration depth measured 
experimentally. We find that X e g can be either smaller or larger than the average penetration 
depth, depending on the correlation length of A. More importantly, the variance and correla- 
tion length of A will likely change with temperature, endowing the experimentally measured 
penetration depth with temperature dependence that is unrelated to the superconducting 
order parameter. 

This work shows that there can be a disorder-induced change of the penetration depth 
that is not caused by the structure of the superconducting gap in momentum space. Rather, 
it reflects the real space variations of the order parameter. An interpretation that assumed 
a spatially uniform penetration depth would infer a larger modulation of A(fc) than truly 
exists. Because A(fc) is the starting point for investigations of the mechanism of the su- 
perconductor, this omission could lead us astray when we seek to determine the underlying 
mechanism. 

How significant is the effect of disorder-induced change in the penetration depth, given 
that A cff / (A) approaches 1 over a large segment of the R(0)-l plane? Modern measurements 
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can routinely resolve sub-nanometer changes in the penetration depth;-^^ in cuprates and 
pnictides the penetration depth is approximately 200 nm, and a 1-nm change in A yields ^ 
of 0.5% - making even small changes in A e fi/ (A) potentially significant. 

Two issues are worth emphasizing. First, we have made no assumption about the dis- 
tribution of A(r), i.e., whether it is normally distributed or follows a different probability 
distribution. However, the calculation presented here only extends to second order, and 
any non-normality only enters at third order and above. Second, A e g has a complicated de- 
pendence on the correlation function R(z), and we know neither its functional form nor its 
temperature dependence. Hence we cannot make any tidy prediction for the low-temperature 
behavior of A(T); there is no power-law to be had. 

Even without perfect knowledge of R(z), it may be possible learn more about A c fr by 
taking advantage of the general constraints that apply to all correlation functions.— 1 ^ In 
particular, the strong similarities between the two cases presented here (Figs. [2] and Hj) lead 
us to expect qualitatively similar behavior in A e fr for most possible correlation functions. 

To make a stronger statement about A(T), we need to determine the local supercon- 
ducting properties of a given chemically doped and intrinsically disordered material, which 
naturally depends on the microscopic details of the superconducting mechanism. Although 
it should be possible to extract a local penetration depth or superfluid density from nu- 
merical methods such as solving the Bogoliubov-de Gennes equations on a lattice, to the 
best of our knowledge this has never been attempted. Several groups have calculated the 
disorder-averaged superfluid stiffness using this approach, for both s-wave^ and d-wave 2 ^ 1 ^ 
models. The full temperature dependence of the disorder-averaged superfluid density can 
also be calculated, 29 but is incomplete, for we have shown that the real space inhomogeneity 
of the superconducting state also contributes to the temperature dependence. 

The larger message is that some measured properties of disordered superconductors will 
not be determined by their disorder averages alone; inhomogeneities can affect the measured 
properties in an experiment-dependent manner. For example, the heat capacity will be 
given by the disorder average because it is additive, but we have seen that the penetration 
depth is non-trivially affected by the disorder. Nonetheless, these two experiments are both 
traditionally interpreted as measuring the same thing - the magnitude of the single-particle 
gap, A(fe). 

These results give a specific example of the potential impact of spatial variation on mea- 
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surements of the penetration depth. With a full consideration of the impact of spatial vari- 
ation on different measured quantities, as well as a complete understanding of how random 
chemical doping gives rise to a non-uniform superconducting state, we will be able to inte- 
grate a complete account of the effects of disorder into our understanding of unconventional 
superconductivity. 
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° (A) (A) 

Position, r 

FIG. 1. Sample realizations of a random penetration depth reveal the influence of the correlation 
length, I, and variance, i?(0). The variance, R(0) = (A 2 ) / (A) 2 — 1, controls the width the penetration 
depth distribution, and the correlation length establishes the characteristic length over which A(r) 
changes. 
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FIG. 2. The effective penetration depth is a strong function of the parameters that characterize the 
distribution of local penetration depths. Here we show the value of A e g- as the correlation length, 
I, and variance, -R(O), run across three orders of magnitude. This figure considers the case of 
squared exponential correlations in the penetration depth; a different case is shown in Fig. HI The 
most important features of this color plot are the range of A e g/(A) and the appearance of values 
both above and below 1. The calculation is valid when R(0) -C 1, but we show the region with 
i?(0) > 0.1 to emphasize the trends seen. Any temperature dependence in I or R(0) will contribute 
to A(T). This temperature dependence is not accounted for by the superconducting gap. 
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FIG. 3. The screening can either be enhanced (A c g < (A)) or suppressed (X c g > (A)), depending 
on the correlation length. The curves for I = 0.1(A) and I = 0.01(A) overlap. 
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FIG. 4. The effective penetration depth for Matern correlations when v = 2 (shown here) has 
strong similarities to Fig. [21 which represents the limiting case where v — > oo. These similarities 
imply that the smoothness of the random penetration depth does not strongly affect X e g . 
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